clear;clc;

modelv=load('BP_cut_model_gas_chimney.mat');

v=(modelv.velp);
Q=(modelv.Qp);

[zz,xx]=size(v);
dx=8;
z_tick=[0:20:zz];z_label=dx*z_tick;
x_tick=[0:50:xx];x_label=dx*x_tick;
Font=14;

figure;
set(gcf,'units','normalized','position',[0.1 0.1 1.0 0.5]);
subplot(1,2,1);
imagesc(v);
ylim([0,zz]);xlim([0,xx]);
set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('x(m)','FontSize',Font);
ylabel('Depth(m)','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'XTick', x_tick);
set(gca,'XTickLabel',x_label,'FontSize',Font);
set(gca, 'YTick', z_tick,'FontSize',Font);
set(gca,'YTickLabel',z_label);

subplot(1,2,2);
imagesc(Q);
ylim([0,zz]);xlim([0,xx]);
set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('x(m)','FontSize',Font);
ylabel('Depth(m)','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'XTick', x_tick);
set(gca,'XTickLabel',x_label,'FontSize',Font);
set(gca, 'YTick', z_tick,'FontSize',Font);
set(gca,'YTickLabel',z_label);
